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Abstract. In this paper we consider the Stochastic Matching problem, which is motivated by appli¬ 
cations in kidney exchange and online dating. We are given an undirected graph in which every edge 
is assigned a probability of existence and a positive profit, and each node is assigned a positive integer 
called timeout. We know whether an edge exists or not only after probing it. On this random graph we 
are executing a process, which one-by-one probes the edges and gradually constructs a matching. The 
process is constrained in two ways: once an edge is taken it cannot be removed from the matching, and 
the timeout of node v upper-bounds the number of edges incident to v that can be probed. The goal is 
to maximize the expected profit of the constructed matching. 

For this problem Bansal et al. [4] provided a 3-approximation algorithm for bipartite graphs, and a 
4-approximation for general graphs. In this work we improve the approximation factors to 2.845 and 
3.709, respectively. 

We also consider an online version of the bipartite case, where one side of the partition arrives node by 
node, and each time a node h arrives we have to decide which edges incident to h we want to probe, and 
in which order. Here we present a 4.07-approximation, improving on the 7.92-approximation of Bansal 
et al. [1]. 

The main technical ingredient in our result is a novel way of probing edges according to a random but 
non-uniform permutation. Patching this method with an algorithm that works best for large probability 
edges (plus some additional ideas) leads to our improved approximation factors. 


1 Introduction 

In this paper we consider the Stochastic Matching problem, which is motivated by applications in kidney 
exchange and online dating. Here we are given an undirected graph G = (U, E). Each edge e € A is labeled 
with an (existence) probability pe G (0,1] and a weight (or profit) We > 0, and each node v € V with a 
timeout (or patience) ty € N^. An algorithm for this problem probes edges in a possibly adaptive order. 
Each time an edge is probed, it turns out to be present with probability pe, in which case it is (irrevocably) 
included in the matching under construction and provides a profit We ■ We can probe at most ty edges among 
the set S{u) of edges incident to node u (independently from whether those edges turn out to be present or 
absent). Furthermore, when an edge e is added to the matching, no edge / € ^(e) (i.e., incident on e) can 
be probed in subsequent steps. Our goal is to maximize the expected weight of the constructed matching. 
Bansal et al. [3] provide an LP-based 3-approximation when G is bipartite, and via reduction to the bipartite 
case a 4-approximation for general graphs (see also [3]). 

We also consider the Online Stochastic Matching with Timeouts problem introduced in [3]. Here we are 
given in input a bipartite graph G = (A U H, A x B), where nodes in B are buyer types and nodes in A are 
items that we wish to sell. Like in the offline case, edges are labeled with probabilities and profits, and nodes 
are assigned timeouts. However, in this case timeouts on the item side are assumed to be unbounded. Then 
a second bipartite graph is constructed in an online fashion. Initially this graph consists of A only. At each 
time step one random buyer b of some type b is sampled (possibly with repetitions) from a given probability 
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distribution. The edges between b and A are copies of the corresponding edges in G. The online algorithm 
has to choose at most tb unmatched neighbors of 6, and probe those edges in some order until some edge ab 
turns out to be present (in which case ab is added to the matching and we gain the corresponding profit) or 
all the mentioned edges are probed. This process is repeated n times, and our goal is to maximize the final 
total expected profilij. 

For this problem Bansal et al. [4] present a 7.92-approximation algorithm. In his Ph.D. thesis Li [8] 
claims an improved 4.008-approximation. However, his analysis contains a mistake [9]. By fixing that, he 
still achieves a 5.16-approximation ratio improving over [4]. 


1.1 Our Results 

Our main result is an approximation algorithm for bipartite Stochastic Matching which improves the 3- 
approximation of Bansal et al. [1] (see Section [5|) . 

Theorem 1. There is an expected 2.8Ab-approximation algorithm for Stochastic Matching in bipartite graphs. 

Our algorithm for the bipartite case is similar to the one from [4] , which works as follows. After solving a 
proper LP and rounding the solution via a rounding technique from [7] , Bansal et al. probe edges in uniform 
random order. Then they show that every edge e is probed with probability at least Xe ■ gipmax), where Xe 
is the fractional value of e, Pmax ■= ^Six.f^s{e){Pf} is the largest probability of any edge incident to e (e 
excluded), and g{-) is a decreasing function with ^(1) = 1/3. 

Our idea is to rather consider edges in a carefully chosen non-uniform random order. This way, we are 
able to show (with a slightly simpler analysis) that each edge e is probed with probability Xe • g (pe) ^ ^Xe- 
Observe that we have the same function g{-) as in [3], but depending on pg rather than Pmax- In particular, 
according to our analysis, small probability edges are more likely to be probed than large probability ones 
(for a given value of Xg), regardless of the probabilities of edges incident to e. Though this approach alone 
does not directly imply an improved approximation factor, it is not hard to patch it with a simple greedy 
algorithm that behaves best for large probability edges, and this yields an improved approximation ratio 
altogether. 

We also improve on the 4-approximation for general graphs in [4] . This is achieved by reducing the general 
case to the bipartite one as in prior work, but we also use a refined LP with blossom inequalities in order to 
fully exploit our large/small probability patching technique. 

Theorem 2. There is an expected 3.709-approximation algorithm for Stochastic Matching in general graphs. 

Similar arguments can also be successfully applied to the online case. By applying our idea of non-uniform 
permutation of edges we would get a 5.16-approximation (the same as in [8], after correcting the mentioned 
mistake). However, due to the way edges have to be probed in the online case, we are able to finely control the 
probability that an edge is probed via dumping factors. This allows us to improve the approximation from 
5.16 to 4.16. Our idea is similar in spirit to the one used by Ma m in his neat 2-approximation algorithm 
for correlated non-preemptive stochastic knapsack. Further application of the large/small probability trick 
gives an extra improvement down to 4.07 (see Section [S]). 

Theorem 3. There is an expected A.07-approximation algorithm for Online Stochastic Matching with Time¬ 
outs. 


1.2 Related work 

The Stochastic Matching problem falls under the framework of adaptive stochastic problems presented first 
by Dean et al. [6]. Here the solution is in fact a process, and the optimal one might even require larger than 
polynomial space to be described. 

The Stochastic Matching problem was originally presented by Chen et al. [S] together with applications in 
kidney exchange and online dating. The authors consider the unweighted version of the problem, and prove 
that a greedy algorithm is a 4-approximation. Adamczyk [T] later proved that the same algorithm is in fact 


^ As in [4], we assume that the probability of a buyer type b is an integer multiple of 1/n. 
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a 2-approximation, and this result is tight. The greedy algorithm does not provide a good approximation in 
the weighted case, and all known algorithms for this case are LP-based. Here, Bansal et al. [4] showed a 3- 
approximation for the bipartite case. Adamczyk [2] presented a different analysis of the same algorithm. Via 
a reduction to the bipartite case, Bansal et al. [4] also obtain a 4-approximation algorithm for general graphs. 
The same approximation factor is obtained by Adamczyk et al. [3] using iterative randomized rounding. 


2 Stochastic Matching 


2.1 Bipartite graphs 

Let us denote by OPT the optimum probing strategy, and let E [OPT] denote its expected outcome. Consider 
the following LP: 


max WePgXg 

(LP-BIP) 


eGE 



S.t. ^ PgXg < 1, 

Vu G V; 

(1) 

eGS{u) 



^ ^ 1 

Vu G V; 

(2) 




0 < Xe < 1, 

Ve G E. 

(3) 

The proof of the following Lemma is alreadv quite standard |3I4I6| 

— just note that Xg 

= P [OPT probes e] 


is a feasible solution of LP-BIP. 


Lemma 1. Let LPup he the optimal value of LP-BIP. It holds that LPup > E [OPT], 

Our approach is similar to the one of Bansal et al. [5] (see also Algorithm [1] in the figure). We solve 
LP-BIP: let x = {xe)eGE be the optimal fractional solution. Then we apply to x the rounding procedure by 
Gandhi et al. [7], which we shall call just GKPS. Let E be the set of rounded edges, and let Xe = 1 if e G A 
and Xe = 0 otherwise. GKPS guarantees the following properties of the rounded solution: 

1. (Marginal distribution) For any e G £1, P [xe = 1] = Xg. 

2. (Degree preservation) For any u G V, J2e&s{v) ^ TSeeiCD) < U. 

3. (Negative correlation) For any v € V, any subset S C 6{v) of edges incident to v, and any b G {0,1}, it 

holds that P [Ae6s(xe = &)] < rieeS ^ • 


Our algorithm sorts the edges in E according to a random permutation and probes each edge e G E according 
to that order, but provided that the endpoints of e are not matched already. It is important to notice that, 
by the degree preservation property, in E there are at most ty edges incident to each node v. Hence, the 
timeout constraint of v is respected even if the algorithm probes all the edges in S{u) n E. 

Our algorithm differs from [3] and subsequent work in the way edges are randomly ordered. Prior work 
exploits a random uniform order on E. We rather use the following, more complex strategy. For each e G E we 


draw a random variable Yg distributed on the interval 


0,^ln 


I-Pe 


according to the following cumulative 


distribution: 'P\Yg<y\ = — (1 — e p®*') . Observe that the density function of Yg in this interval is e 


(and zero otherwise). Edges of E are sorted in increasing order of the Yg’s, and they are probed according 
to that order. We next let Y = {Yg)^^p,. 

Define d{v) := S (v) O E. We say that an edge e G E is safe if, at the time we consider e for probing, no 
other edge / G S{e) is already taken into the matching. Note that the algorithm can probe e only in that 
case, and if we do probe e, it is added to the matching with probability pg. 

The main ingredient of our analysis is the following lower-bound on the probability that an arbitrary 
edge e is safe. 


Lemma 2. For every edge e it holds that P 


e is sa/e| e G E 


> g [pg), where 


9 (P) ■= 


1 


— exp 


(2 -I- p) - In 
P 



2+p 
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Algorithm 1 Approximation algorithm for bipartite Stochastic Matching. 

1. Let (®e)gg£ be the solution to LP-BIP. 

2. Round the solution with GKPS; let {xs)s£e be the rounded 0-1 solution, and E = {e £ E\xe = 1}. 

3. For every e € E, sample a random variable Yg. distributed as P \Yg. < y] = ^ ■ 

4. For every e € F in increasing order of Fei 

(a) If no edge / € S{e) := 5{e) D E is yet taken, then probe edge e 


Proof. In the worst case every edge / G (5(e) that is before e in the ordering can be probed, and each of these 
probes has to fail for e to be safe. Thus 


e is safe| eGE 

— ®".E\e,V 

n o-pf) 

eGE 



j€S{ey.Yf<Y, 



Now we take expectation on Y only, and using the fact that the variables Yf are independent, we can write 
the latter expectation as 






In , 


Y[ (P [Yf <y]{l- Pf) +F[Yf>y])\ 


e£ E 


(4) 


Observe that P [Yf < y] (1 — pf) -I- P [Y/ > y] = 1 — p/P [Yf < y] . When y > ^ In , then P [Yf < y] = 
1, and moreover, ^(1 ~ is an increasing function of y. Thus we can upper-bound P [Y/ < y] by 

^(1 — for any y G [0, c»], and obtain that 1 —p/P [Yf < y] > 1 —p/^(l — . Thus (|4]) 

can be lower bounded by 


^E\e 


—^E\e 


f Pe 1 —Pe ^ 

I g ^/G< 5 (e) Pf y P^ ^diy 

Jo 

_ l _ (l- e-(S/eS(e) P/+Pe)it In 

Y.f^S(e)Pf+Pe^ ' 


e£ E 


eG E 


From the negative correlation and marginal distribution properties we know that 


Xf\e G E 


< 


®£)\e [*/] = f ^ He), 


and therefore E^^g 


E/ei(e) Pf 


eG E 


^ J 2 fe 5 {e)Pf^f ^ 2, where the 


last inequality follows from the LP constraints. Consider function f{x) := ~ ® i^+Po) p,. i" i-pe ^. This 

function is decreasing and convex. From Jensen’s inequality we know that E [f{x)] > /(E [a:]). Thus 


E 


E\e 


/1 E p/ 

JeHe) 


e G E 


> / I E^\e 
1 


E Pf 

.f£S{e) 


e G E 


- ^ ® 1-Pe 'j = g[pg.). 

2+Pe\ J 


□ 


From Lemma [5] and the marginal distribution property, the expected contribution of edge e to the profit 
of the solution is 


WePe • J 


eGE 


e is safe] eGE > WePeXe ■ g{Pe) > WePeXe ' 5(1) = ^WePeXe- 


Therefore, our analysis implies a 3 approximation, matching the result in [4]. However, by playing with the 
probabilities appropriately we can do better. 
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Patching with Greedy. We next describe an improved approximation algorithm, based on the patching of the 
above algorithm with a simple greedy one. Let 5 € (0,1) be a parameter to be fixed later. We define Eiarge 
as the (large) edges with pe ^ <5) and let be the remaining (small) edges. Recall that LPi,ip denotes 

the optimal value of LP-BIP. Let also LPiarge and LP^^aii be the fraction of LP^ip due to large and small 
edges, respectively; i.e., LPiarge = Sees,arse '^ePeXe and LP^maii = LPb^p - LPiarge- Define 7 G [0,1] such 
that ^LPbip = LPiarge- By refining the above analysis, we obtain the following result. 

Lemma 3. Algorithm\^ has expected approximation ratio ^7 + g(S) (1 — 7 ). 

Proof. The expected profit of the algorithm is at least: 


WePeXe ' g(Pe) > WePeXe ' g(^) P WePeXe ' g(S) 


= ^LPiarge + g(S)LPsmall = ^7 + g(^) (1 “ 7 )^ 


LPi 


bip • 


□ 


Consider the following greedy algorithm. Compute a maximum weight matching Mgrd in G with respect to 
edge weights WePe, and probe the edges of Mgrd in any order. Note that the timeout constraints are satisfied 
since we probe at most one edge incident to each node (and timeouts are strictly positive by definition and 

w.l.o.g.). 

Lemma 4. The greedy algorithm has expected approximation ratio 5^. 


Proof. It is sufficient to show that the expected profit of the obtained solution is at least 6 - LPiarge- Let 
X = (xe)eGE be the optimal solution to LP-BIP. Consider the solution x' = (x'^fe^E that is obtained from x 
by setting to zero all the variables corresponding to edges in Esmall, and by multiplying all the remaining 
variables by 5. Since Pe > S for all e G Eiarge, x' is a feasible fractional solution to the following matching 
LP: 


max WePeZe (LP-MATCH) 

eeE 

s.t. Ze < 1, VuGR; 

e£S{u) 

0<Ze<l, Ve G R. (5) 

The value of x' in the above LP is d - LPiarge by construction. Let LPmatch be the optimal profit of LP- 
MATCH. Then LPmatch > 5 - LPiarge- Given that the graph is bipartite, LP-MATCH defines the matching 
polyhedron, and we can find an integral optimal solution to it. But such a solution is exactly a maximum 
weight matching according to weights WePe, i-e- X^eeM ^ '^ePe = LPmatch- The claim follows since the 
expected profit of the greedy algorithm is precisely the weight of Mgrd- Cl 

The overall algorithm, for a given d, simply computes the value of 7 , and runs the greedy algorithm if 
71 ^ > (^7 + g(d) (1 — 7 )), and Algorithm [T] otherwis^. 

The approximation factor is given by max{^ + (1 — 7 ) 5 (( 5 ), 7 ^}, and the worst case is achieved when the 

two quantities are equal, i.e., for 7 = 1 , yielding an approximation ratio of ' Maximizing 

(numerically) the latter function in d gives d = 0.6022, and the final 2.845-approximation ratio claimed in 
Theorem [TJ 


2.2 General graphs 

For general graphs, we consider the linear program LP-GEN which is obtained from LP-BIP by adding the 
following blossom inequalities: 

yy PeXe < ^ VW C R |W| odd. ( 6 ) 

e^E(W) 

Note that we cannot run both algorithms, and take the best solution. 


5 






6 


Here E{W) is the subset of edges with both endpoints in W. We remark that, using standard tools from 
matching theory, we can solve LP-GEN in polynomial time despite its exponential number of constraints; 
see the book of Schrijver for details m- Also in this case Xe = P [OPT probes e] is a feasible solution of 
LP-GEN, hence the analogue of Lemma [T] still holds. 

Our Stochastic Matching algorithm for the case of a general graph G = (V, E) works via a reduction to 
the bipartite case. First we solve LP-GEN; let x = (xe)gg^ be the optimal fractional solution. Second we 
randomly split the nodes V into two sets A and B, with Eab being the set of edges between them. On the 
bipartite graph (AU B, Eab) we apply the algorithm for the bipartite case, but using the fractional solution 
{xe)e^EAB ii^cluced by LP-GEN rather than solving LP-BIP. Note that {xe)^^EAB ^ feasible solution to 
LP-BIP for the bipartite graph {A U 5, Eab)- 

The analysis differs only in two points w.r.t. the one for the bipartite case. First, with Eab being 


the subset of edges of Eab that were rounded to 1 , we have now that 


e G Eab 


— P [e G Eab] 


e G Eab 


e G Eab 


= hxe- Second, but for the same reason, using again the negative correlation and 


marginal distribution properties, we have 


E 

1 

M 

e G Eab 





< 


E pf^ 

feS{e) 


I^Eab 


\ E pf^f - 

fesie) 


Repeating the steps of the proof of Lemma [5] and including the above inequality we get the following. 


Lemma 5. For every edge e it holds that P 


e is safe\ e G Eab 


> h{pe), where 


* w TT; (i - “P (- (1+f) ;1" tv;)) ■ 

Since h[pe) > h{l) = i, we directly obtain a 4-approximation which matches the result in [4]. Similarly 
to the bipartite case, we can patch this result with the simple greedy algorithm (which is exactly the same 
in the general graph case). For a given parameter S G [0,1], let us define 7 analogously to the bipartite 
case. Similarly to the proof of Lemma [Sj one obtains that the above algorithm has approximation factor 
^ -|- ^-^h{6). Similarly to the proof of LemmalU the greedy algorithm has approximation ratio jS (here we 
exploit the blossom inequalities that guarantee the integrality of the matching polyhedron). We can conclude 
similarly that in the worst case 7 = 2 S+h{S)-i /2 > yielding an approximation ratio of 2 S+h\&)-i /2 ■ Maximizing 
(numerically) this function over 5 gives, for 5 = 0.5580, the 3.709 approximation ratio claimed in Theorem 

H 


3 Online Stochastic Matching with Timeouts 

Let G = (A U R, A X R) be the input graph, with items A and buyer types R. We use the same notation for 
edge probabilities, edge profits, and timeouts as in Stochastic Matching. Following [4], we can assume w.l.o.g. 
that each buyer type is sampled uniformly with probability 1/u. Gonsider the following linear program: 


max ^ WabPabXab 

(LP-ONL) 

a£A,b£B 


S-t. ^ ^ Pab^ab E I 5 

Va G A 

bGB 


^ ^ Pab^ab ^ I 5 

yb€B 

a£A 


^ ^ ^ab ^ ^67 

yb€B 

a£A 


0 ^ ^ab ^ I 7 

Vafc G E. 
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The above LP models a bipartite Stochastic Matching instance where one side of the bipartition contains 
exactly one buyer per buyer type. In contrast, in the online case several buyers of the same buyer type (or 
none at all) can arrive, and the optimal strategy can allow many buyers of the same type to probe edges. 
Still, that is not a problem since the following lemma from [J allows us just to look at the graph of buyer 
types and not at the actual realized buyers. 

Lemma 6. Lemmas 9 and 11) Let E [OPT] be the expected profit of the optimal online algorithm for 
the problem. Let LPoni be the optimal value of LP-ONL. It holds that E [OPT] < LPoni- 

We will devise an algorithm whose expected outcome is at least ■ LPoni, and then Theorem |3] follows 
from Lemma ini 

The algorithm. We initially solve LP-ONL and let (xab)abeAxB optimal fractional solution. Then 

buyers arrive. When a buyer of type b is sampled, then 1) if a buyer of the same type b was already sampled 
before we simply discard her, do nothing, and wait for another buyer to arrive, 2) if it is the first buyer of 
type 6, then we execute the following subroutine for buyers. Since we take action only when the first buyer 
of type b comes, we shall denote such a buyer simply by b, as it will not cause any confusion. 


Subroutine for buyers. Let us consider the step of the online algorithm in which the first buyer of type 
b arrived, if any. Let At be the items that are still available when b arrives. Our subroutine will probe a 
subset of at most tb edges ab, a € Ab. Consider the vector {xab)a^Ah - Observe that it satisfies the constraints 
J2a(^AiPabXab < 1 and J2a(^Ai ^ab < 4- Again using GKPS, we round this vector in order to get 
with Xab £ and satisfying the marginal distribution, degree preservation, and negative correlation 

propertied. Let Ab be the set of items a such that Xab = 1- For each ab, a G Ab, we independently draw 
a random variable Yab with distribution: P [Yab < y] = — (1 — exp {—pab ■ y)) for y £ 


0, — In y- 

’ Pa.b 1- 


Let 


y = iyab)o^Ab- 

Next we consider items of Ab in increasing order of Yab. Let aab £ [^, 1] be a dumping factor that we 
will define later. With probability aab we probe edge ab and as usual we stop the process (of probing edges 
incident to b) if ab is present. Otherwise (with probability 1 — aab) we simulate the probe of ab, meaning 
that with probability pab we stop the process anyway — like if edge ab were probed and turned out to be 
present. Note that we do not get any profit from the latter simulation since we do not really probe ab. 


Dumping factors. It remains to define the dumping factors. For a given edge ab, let 

fiab ■= 

Using the inequality X]ae. 4 i, by repeating the analysis from Section [5] we can show that 

I 


n (l-Pa'b) 

aG Ab 

a' eAb:Y^,,,<Ya^b 



Pab > h{pab) = 


I +Pab 


1 - exp ( - (I -I- Pab) — In -— - - 

Pab 1 - Pab 


1 

> -. 

- 2 


Let us assume for the sake of simplicity that we are able to compute Pab exactly. We will show in Appendix 
m how to remove this assumption. We set a^b = defined since Pab £ [1/2, !]• 


Analysis. Let us denote by Ab the event that at least one buyer of type b arrives. The probability that an 
edge ab is probed can be expressed as: 

P [Mb] • P [no b' takes a before b] Mb] • P [b probes a] Mb A a is not yet taken]. 

The probability that b arrives is P [Mb] = 1 — (l — L])” >1 — 4 We shall show first that 

P [b probes a] Mb A a is not yet taken] 

® Actually in this case we have a bipartite graph where one side has only one vertex, and here GKPS reduces to 
Srinivasan’s rounding procedure for level-sets m- 











is exactly hxab: and later we shall show that P[no b' takes a before b\Ab] is at least — j-r —rr- This will 

1+2 11 “ e j 

yield that the probability that ab is probed is at least 

/ 1\ 1 1 e-1 1 

[ ej 1 + 1(1-i) ■ 2^“'“ 3e-l^“'^ 4.16'"“''- 


Consider the probability that some edge a'h appearing before ah in the random order blocks edge ab, 
meaning that ab is not probed because of a'b. Observe that each such a'b is indeed considered for probing in 
the online model, and the probability that a'b blocks ab is therefore Oia'bPa'b + (1 ~ cta'b)Pa'b = Pa'b- We can 
conclude that the probability that ab is not blocked is exactly /3o&- 

Due to the dumping factor aab, the probability that we actually probe edge ab € Ab is exactly aab'fiab = \- 


Recall that P 


a G Ab 


Xab by the marginal distribution property. Altogether 


P [6 probes a| A& A a is not yet taken] = -Xab- 


(7) 


Next let us condition on the event that buyer b arrived, and let us lower bound the probability that ab is 
not blocked on the a’s side in such a step, i.e., that no other buyer has taken a already. The buyers, who are 
first occurrences of their type, arrive uniformly at random. Therefore, we can analyze the process of their 
arrivals as if it was constructed by the following procedure: every buyer b' is given an independent random 
variable Yb' distributed exponentially on [0, oo], i.e., P [Ybi < y] = 1 — buyers arrive in increasing order of 
their variables Yb'. Once buyer b' arrives, it probes edge ab' with probability (exactly) o-awPab'Xab' = \xab' 
— these probabilities are independent among different buyers. Thus, conditioning on the fact that b arrives, 
we obtain the following expression for the probability that a is safe at the moment when b arrives: 


[no b' takes a before b\ Ab] 

> E {1-¥[Ab'\Ab]'P[b'probes ab'\Ab']pab') 

b'eB\b-.Y^,<Yb 


Ab 


poo 

/ TT [l - F [Ab'\ Ab] ■'P [Yb' <y\Ab']-P[b'probes ab'\Ab']pab')e~ydy. 


b’eB\b 


Now let us upper-bound each of the probability factors in the above product. First of all P[A6'|A6] = 
1 — (l — i)” ^ < 1 — Second, P [Yb' < y\ Ab'] = 1 — just by definitioil^. Third, from ([7]) we have that 
P [&' probes ab'\ Ab'] = 

Thus the above integral can be lower bounded by 



- e ^) ■ ^Xab' -Pab'^ e ^dy 
^Xab' -Pab' ■ yj e“^dy 


> 


1 

1 + (l — ^) 5 (^b'eB\bPab' • Xab''^ 

1 _ 2e 


Above in the first inequality we used the fact that 1 — c(l — e ^) > e for c S [0,1] and any y € M.: here 
c = (l — i) ^Xab' 'Pab'- In the first equality we used e~°'^dx = A In the last inequality we used the LP 
constraint J2b'eB\bPab' ■ Xab' < 1- 

^ The Ab' event in the condition simply indicates that Yb' was drawn. 
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Altogether, as anticipated earlier, 


P [ai, is probed] > (l - i) ^ ■ ^ > ■ x.s. 

In Appendix [B] we will show how to compute the dumping factors so that the above probability is 
for an arbitrarily small constant e > 0. In particular, by choosing a small enough e the factor 4.16 is still 
guaranteed. 

We can again use the approach with big and small probabilities, thus reducing the approximation factor 
to 4.07. The details are given in Appendix lAl Theorem |3] follows. 
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A Combination with Greedy in the Online Case 

Recall that h{p) = ^1 — exp (1 + p) Mn We are again applying the big/small probabilities 

trick, so let 5 £ (0,1) be a parameter to be fixed later. Consider back again the subroutine for buyers. 
Previously we have used dumping factors Uab = — recall — flab > h {pab)- 

This time we define Uab = for such that pab < S, and aab = otherwise. We again assume 

here that we can calculate Pab (see Appendix I bI). Define Eiarge = {ab £ E \pab > and Esmaii = E\Eiarge, 
and let LPiarge = 7 • LPoni- Therefore, for edge ah the probability that ab is probed when b scans items is 
exactly h (d) for ab £ Es^naii and i for ab £ Eiarge ■ Now by repeating the steps in the proof of Section [31 we 
obtain that the probability that ab is not blocked on a’s side is at least 

1 ^ 1 

1 T (1 ~ e) i^hb' ^B\bPab’ ' O^ab' Pab' ‘ Xab'^ ^ ~ e) ^ i^hb' ^B\bPab' ' Xab'^ 

1 

-1 + { 1 -I)hisy 
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since Uab' ’Pab' = h (S) for small edges and aab' 'Pab' = ^ < h{S) for large edges. Therefore, the approximation 
ratio of such an algorithm is at least 

(i-l) (, _V2_^ \ 

An alternative algorithm simply computes a maximum weight matching w.r.t. weights PeWe in the graph 
corresponding to LP-ONL, and upon arrival of the first copy of a buyer type b probes only the edge incident 
to b in the matching (if any). By the same argument as in the offline case, this matching has weight at least 
'y ■ S ■ LPoni , and every buyer type is sampled with probability at least 1 — . So the approximation ratio of 

the greedy algorithm is at least (l — i) yj. 

For a fixed 5, depending on the value of 7 (that we can compute offline) we can run the algorithm with 
best approximation ratio according to the above analysis. Thus the overall approximation ratio is 

Optimizing over 6 gives 5 = 0.525 and a final approximation factor strictly less than 4.07. 

B Computing Dumping Factors 

Recall that we assumed the knowledge of quantities Pab-, which are needed to define the dumping factors 
Uab- Though we are not able to compute the first quantities exactly in polynomial time, we can efficiently 
estimate them and this is sufficient for our goals. Let us focus on a given edge ab. Recall that 


n (l-Pa'b) 

aG Ab 




> - (1 - exp ( - (1 + Pah) — In—-) ) =h[pab) ■ 

1 + Pah V V Pah 1 - Pah J ) 


Let us simulate the subroutine for buyers N times without the dumping factors —- in a simulation we 
run GKPS, we sample the Y variables, but we simulate probes of edges, and we never really probe any edge. 
We shall set lY later. Let be 0-1 indicator random variables of whether a was safe or not in 

each simulation. Note that E [S'®] = PabXab S [h (pab) Xab, Xab]- 

Suppose that Xab > where n is the number of buyers. The expression Sab = should be a good 

estimation of Pab-Xab, i-e.. Sab € [PabXab (1 - e), PabXab (1 + e)] with probability 1-;^. Set ^ In {2rPZ) 
for Z = 3^ + 1. 

Applying Chernoff’s bound P [|A - E [A] | > eE [A]] < 26“^with A = -5, one obtains: 


N 


^ Sj ^ [(1 - e) PabXab • A, (1 -I- e) PabXab ' N] 




<2 exp ——Pah^ah * < 2 exp ( — — < 2 exp —— • N 


(P Xab 


3 2 


6 n 


1 1 


From the union-bound, with probability at least 1 — A we have that Sab S [PabXab (1 — e)) PabXab (1 + e)] 
for every edge ab such that Xab > 

Now let us assume this happened, i.e., we have good estimates. We set aah = max{i, min{i , 1}} which 


belongs to 


1 1 1 1 
2 /3ai,(H-e)’ 2 /3„b(l-e) 


but only for edges ab such that Xab > -■ For edges ab such that Xab < — 
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we just put Uab = 1 (so we do not dump such edges actually). Two elements of the proof were depending on 

Xgb _Xab_ 

2(l+e)’ 2(l-e) 


the dumping factors. First, now the probability that edge is taken is aabPabXab G i ■ Second, 


recall the probability of an edge ab not to be blocked: 

1 


We have that 


< 


1 T (l g) {^2ib'^B\bPo^b’ ‘ Oiab'Pab' ' Xab’^ 


^ ^ Pab' ‘ ^ab' i^ab' * ^ab' 
b'GB\b 

^ ^ Pab' * ^ab'Pab' * ^ab' T ^ ( Pab' ‘ ^ab'Pab' ' ^ab' 

b'^B\b-.x^y>j^ b'£B\b-.x^y<j; 

1 


( 8 ) 


E 


Pab' ■ 


2(l-e) 


^ab' 


E 


^ab' 


b'eB\b-.x^^,>-^ 




So the probability that a is not blocked is at least — i)(i+o( )) • final probability that edge ab is 

probed is at least 


1 

1 - - 


^ab 


^ab 


e — 1 


e;2(l + e) l+(l-i) (i + 0(e)) l + £ 2e + (e - 1) (1 + O (e)) 


— ^ab 


e — 1 


> 


1 


3e-l + 0(e) 4.16 


^ab’ 


In the last inequality above we assumed e to be small enough. 

With probability at most we did not obtain good estimates of the dumping factors. Still we have that 
ctab G l] I therefore UabPab G [ji l] • this case quantity ([5]) can be just lower-bounded by 

and the probability that edge ab is probed in the subroutine for buyers is at least Thus the probability 
that edge ab is probed during the algorithm is at least (l — |) ^ ^ > Yjygja^ah- The 

total expected outcome of the algorithm is therefore, for sufficiently small e, at least 


LP, 


onl 


1 

1 - ^ 


e — 1 


Z J 2,6-1 +0{e) Z A 2e - 1 


11 e — 1 \ ^-3^+1 1 

> —rLPonl- 


4.16 


The above approach can be combined with the small/big probability trick from Appendix!^ By choosing 
e small enough the approximation ratio is 4.07 as claimed. 























